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1.0  INTRODUCTION 

CNO  Operational  Requirement  XCC20  of  23  June  1975  establishes  the 
need  for  low  probability  of  intercept  communications  for  installation 
aboard  submarines  and  aircraft.  In  addition  to  covertness,  a communica- 
tion system  is  required  to  be  independent  of  the  submerged  vehicle's 
speed  and  depth.  Generally,  submarines  are  antenna-limited  in  their 
communications.  While  several  radio  frequency  systems  are  presently 
operational  or  undergoing  OPEVAL/TECHEVAL,  their  use  severely 
handicaps  the  submarine's  speed  and  depth  in  direct  support  of  surface 
units.  While  it  is  recognized  that  no  single  system  can  meet  all  Fleet 
requirements,  only  a submar ine-to-aircraf t optical  communication  system 
can  achieve  absolute  covertness  and  speed  independence.  However,  such 
a system  is  generally  assumed  to  be  only  effectual  in  cloud-free 
environments;  a limitation  which  greatly  impairs  the  overall  system 
capability.  Vet  the  fact  that  sunlight  gets  through  optically  dense 
clouds  on  overcast  days  suggests  that  a sufficient  amount  of  pulsed 
laser  radiation  might  be  available  for  demodulation/decoding  after  trans- 
versing  the  cloud  if  one  made  the  f ield-of-view  (FOV)  of  the  receiver 
wide  enough  to  collect  the  scattered  irradiance.  The  effects  of  multiple 
scattering  by  dense  atmospheric/marine  clouds  on  pulsed  radiation  trans- 
lates into  a pulse  stretching  and  spatial  broadening  of  the  initial 
signal.  Thus  the  performance  of  any  gated  and  continuous  image  viewing 
systems  and  communications  systems  using  pulse-coding  techniques  will  be 
degraded  in  terms  of  limited  message  data  rates  and  reduced  signal-to- 
noise  ratios. 

The  Navy  presently  does  not  have  sufficient  information  about  the 
optical  propagation  properties  of  aerosols  in  the  marine  boundary  layer 
to  provide  designers  and  users  of  electrooptics  (EO)  systems  with  parameters 
required  to  determine  optimum  system  design  and  hence  operational  system 
limitations . 

Measurements  of  the  direct  beam  extinction  coefficients  (absorption 
plus  single  scatter)  of  optical  frequencies  through  various  continental 
type  clouds,  fogs,  hazes  and  rain  have  been  widely  reported  in  the  open 
literature.  Recently,  Bucher  and  Lerner  [22]  reported  measurements 
of  multiple  scattering  effects  on  30  nsec  pulses  from  a Q-switched  ruby 
laser  propagating  through  orographic  clouds.  The  pulses  received  with 
varying  FOV  (1°  to  4°)  over  a 6.4  km  path  varied  between  1 and  10  usee 
in  duration  when  there  was  a cloud  in  the  propagation  path.  Quantitative 
measurements  of  the  peak  power  losses  as  a function  of  receiver  FOV  were 
not  addressed  and  detailed  descriptions  of  the  clouds  (aerosol  and 
hydrometer  size  distributions  and  their  spatial  variations)  were  not 
-.available.  ..  f 

\ ' 

"“it  is  the  purpose  of  this  test  plan  to  define  an  experiment  to  measure 
the  peak  power  losses  and  temporal  characteristics  as  a function  of  FOV  of 
optical  pulses  propagating  through  clouds  (marine  fogs) . In  addition  to  the 
propagation  measurements  the  meteorological  parameters  describing  the 
propagation  path  will  be  measured  and  used  as  inputs  to  computer  codes  to 
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assess  the  measured  results  and  for  propagation  prediction  purposes.  Of 
primary  interest  here  are  those  frequencies  in  the  blue-green  portion  of  the 
spectrum  which  after  propagating  through  a cloud,  can  penetrate  sea  water  and 
be  of  use  as  a submar ine-to-aircraft  satellite  communications  system.. 

Recently,  the  Naval  Electronics  Laboratory  Center  in  collaboration  with 
the  Naval  Research  Laboratory  initiated  an  experimental  program  to  investi- 
gate the  propagation  of  infrared  (10.6  w)  radiation  through  optically  thick 
clouds  to  assess  the  usefulness  of  heterodyne  techniques  for  surface-to- 
satellite  communications  [2].  Unfortunately,  this  program  was  terminated 
just  prior  to  data  acquisition  due  to  a recall  of  funds.  Subsequent 
measurements  were  made  over  an  11-mile  link  transversing  San  Diego  Bay 
and  the  results  have  been  reported. [2]  The  instrumentation  and  experi- 
mental expertise  developed  for  this  program  and  other  related  laser  radar 
projects  are  currently  available  for  conducting  the  measurements  proposed 
here.  In  addition,  a unique  multiremote  sensing  facility  has  been 
developed  by  NELC  and  is  available  for  conducting  the  propagation  parameter 
measurements.  In  the  following  sections  the  details  of  the  test  plan 
(site  selection,  facilities  and  instrumentation  description,  variables 
to  be  measured  and  schedule,  theoretical  background  and  results,  and 
computer  software  capabilities)  are  presented. 

2.0  THEORETICAL  BACKGROUND  AND  PREDICTIONS 

2.1  Introduction 

The  existence  of  clouds  and  fog  in  many  regions  of  the  earth  presents 
a formidable  problem  to  the  designer  of  an  optical  communication  system 
whose  transmission  channel  is  the  atmosphere.  These  entities  inhibit 
system  performance  by  inducing  beam  spread,  dispersion  in  angle-of- 
arrival,  degradation  of  spatial  coherence,  and  dispersion  in  time  and 
frequency  of  the  signal  modulating  the  optical  beam.  Thus  one  is 
forced  in  many  applications  to  accept  receiver  complexity  and  diminished 
data  rates  in  return  for  nearly  constant  channel  availability.  The 
purpose  of  this  section  is  to  review  the  various  theoretical  formula- 
tions currently  used  to  characterize  optical  propagation  in  single  and 
multiple  scattering  media. 

2.2  Mie  Scattering  and  the  Radiative  Transfer  Equations 


The  following  discussion  is  a summary  of  the  Radiative  Transfer 
and  Mie  Scattering  Theories  taken  from  a review  paper  by  Lerner  and 
Holland  [3]1. 

2.2.1  Radiative  Transfer:  Preliminary  Considerations 

The  basis  of  classical  transfer  theory  is  the  Lambert-Bouger  law, 
which  states  that  the  change  in  intensity  of  radiation  due  to  an  inter- 
action with  matter  depends  directly  on  both  the  incident  intensity  and 
the  amount  of  matter.  This  linearity  holds  so  long  as  the  matter 
remains  in  the  same  physical  state  and/or  if  we  consider  quantities 
at  a single  frequency  or  wavelength.  Thus  a beam  or  pencil  of  radiation 
is  attenuated  by  its  interaction  with  the  medium  through  which  it  passes. 
If  the  initial  intensity  is  I,  then  after  traveling  a distance  ds  along 
the  direction  of  propagation,  the  change  in  intensity  will  be2 
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density  (. g/crn 3)  of  the  attenuating  medium 
k?  i mass  extinction  coefficient  (cm^Ao) 
f = volume  extinction  coefficient  (cm-^) 


In  what  follows,  only  the  variations  of  intensity  with  resDect  to 
direction  and  spatial  position  will  be  indicated.  The  dependence  on 
frequent'*,  tire,  e:'e-  necessarv  parameters  will  nevertheless  be 

understood  unless  otherwise  stated. 

The  theorv  also  assumes  that  emission  bv  an  element  of  mass  or 
volume  is  linear  in  both  intensity  and  the  amount  of  matter  present. 
The  net  increase  in  intensity  in  a distance  ds  along  the  direction  of 
propagation  is  j(slds.  The  equation  of  radiative  transfer  for  a given 
path  through  a medium  can  he  thus  expressed  as 

dl(s) 

= "6ext(s)Us)  + j(s)  (1) 

Formal  integration  of  this  one-dimensional  transfer  eauation  of  the 
path  s = 0 to  s = T gives 
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j (s  *)e-(T  _ T' } ds ' 


in  which  the  optical  thickness  along  the  path  is  t 
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flpvr  (s”)d s". 


This  formal  solution  to  (1)  shows  that  the  resultant  intensitv  is  the 
sur  of  two  distinct  terms;  the  incident  intensity  attenuated  bv  extinc- 
tior.  along  the  path  (0,T)  and  a component  due  to  emission  at  points 
along  the  path,  attenuated  along  that  path  (s',T).  If  we  neglect  the 
enission  term,  we  have  the  usual  form  of  the  Lambert -Bouger  law.  Thus 
we  can  define  a "meteorological  range"  as  the  distance  required  to  reduce 
the  incident  intensitv  to  ?’(■  of  its  original  value:  that  is, 
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This  relation  is  widely  used;  however,  it  is  deceptive  in  that  the  above 
describes  nothing  about  the  deposition  of  the  light  removed  from  the  beam. 
That  information  is  contained  in  the  emission  term.  In  general,  the 
emission  term  consists  of  two  distinct  parts:  (1)  that  due  to  the 
scattering  of  other  radiation  into  the  direction  of  interest,  and  (2) 
that  due  to  true  emission,  e.g.,  thermal.  Henceforth  we  will  ignore 
this  latter  term,  and  consider  only  emission  due  to  scattering. 

Equations  (1)  or  (2)  are  known  as  the  "equation  of  radiative  transfer”. 
They  serve  to  demonstrate  the  dependence  of  the  classical  theory  on  the 
scattering  centers,  or  on  the  volume  elements  containing  such  scatters. 
Naturally,  these  volume  elements  should  be  small  enough  so  they  may  be 
considered  point  sources  of  radiation;  and  the  number  of  scatterers 
present  ir.  each  volume  not  so  large  that  multiple  scattering  effects  are 
observed.  Under  these  conditions,  one  is  able  to  perceive  a significant 
part  of  the  physics  behind  the  scattering  phenomenon.  It  is  possible 
to  analyze  multiple  scattering  problems  with  this  formation,  e.g., 
Dell-Imagine  [A],  but  the  mathematical  and  computational  complexities 
increase  significantly  unless  a modified  computing  technique,  e.g., 

Dave  and  Gazdag  fllj  or  a matrix  approach,  e.g.,  Plass,  et  al . [5],  is 
utilizer,  "or  this  development,  we  will  limit  ourselves  to  the  single 
scattering  problem. 

In  general , a medium  may  exhibit  both  scattering  and  absorption. 

Some  of  Che  energy  removed  from  the  beam  will  be  scattered  over  all 
possible  directions;  the  rest  will  be  absorbed  and  lost  to  the  radiation 
field.  ■"Vi?  we  have 


kext  ksca  + kabs 


where  the  subscripts  denote  scattering  and  absorption.  An  important 
auxiliarv  parameter  is  the  "albedo  for  single  scattering"  defined  to  be 
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The  mass  scattering  coefficient  specifies  the  rate  at  which  energy 
is  scattered  from  the  incident  beam  confined  to  solid  angle  du  bv  an 
element  of  mass  dm 


dE 

It 


sea  I dm  du> 


k I dm  do 
0 ext 


(3) 


The  distribution  of  the  scattering  energy  is  specified  bv  the  "scalar 
phase  function”  P(V)  where  V is  the  angle  between  the  direction  of 
incidence  and  the  scattered  direction.  Thus,  the  rate  at  which  energy 
is  scattered  into  a solid  angle  dfi  in  a direction  making  an  angle  9 with 
the  direction  of  incidence  of  a beam  of  intensity  I confined  to  a solid 
angle  d>o  from  an  element  of  mass  dm  is 
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dt 
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The  rate  at  which  energy  is  lost  from  the  beam  due  to  scattering  in  all 
directions  is 


k I dmdu  / PCY)  ~ 

SC  a / HIT 

n 
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(5) 


It  is  apparent  that  equations  (3)  and  (5)  agree  only  if  the  scalar 
phase  function  is  normalized  to  one  over  the  unit  sphere,  i.e.. 


POO 


d:"  = 1.0. 


These  quantities  completely  specify  the  single  scattering  properties  of 
[7fSand°?8]  mediunu  For  a more  ^tailed  treatment,  see  references  [6], 


3.2.:  Single  Particle  Scattering  Properties 

«mhpr!le  ^aCt*rin«  a"d  absorption  properties  of  an  isotropic,  homogeneous 
sphere,  illuminated  by  a plane  monochromatic  electromagnetic  wave,  can  be 
described  completely  bv  the  well-known  Mie  theory  [9-11).  We  will  onlv 
give  a brief  description  here.  y 
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Essentially,  given  a diameter  D of  the  spherical  particle,  its 
conplex  refractive  Index  relative  to  the  surrounding  medium  ra  and  the 
wavelength  A of  the  incident  plane  wave,  the  theory  yields  the  dimension- 
less complex  scattering  amplitudes  for  any  direction  8 from  the  direc- 
tion of  incidence.  The  expressions  have  their  simplest  form  when  the 
electric  vector  of  both  incident  and  scattered  waves  are  resolved  into 
components  (1)  perpendicular,  and  (2)  parallel  to  the  plane  of  scattering; 
that  is,  the  plane  defined  by  the  directions  of  incidence  and  scattering. 
cigure  1 depicts  the  geometry  of  the  problem  for  a particular  incident 
polarization.  The  scattering  amplitudes  are  then 


s (a.x.e)  * y 

1 n«l 

2n+l 

[an(m,x)*n(e)  + bn(a,x)£n(e) ] 

n(irH) 

00 

s?(a,x,e)  « y 

2n+l 

fan(fi,x)^n(8)  + bn(ffl,x)irn(6)] . 

n(n+l) 

n=l 


The  quantities  an(m,x)  and  bn(5,x)  are  the  Mie  scattering  coefficients, 
which  are  certain  3icatti~3e=sel  functions  of  real  and  complex  argument. 
"„(?)  and  £ (6)  are  spherical  harmonics  expressable  as  Legendre  poly- 
nomials. The  quantity  x is  called  the  normalized  Mie-size  parameter  and 
is  given  by 


The  components  of  the  scattered  intensity  resolved  in  the  directions 
perpendicular  and  parallel  to  the  scattering  plane  are 

il(m,x,°)  = !Sj(m,x,8>  j 

i2(m,x,9)  = )S2(ffl,x,9)  !2 

The  intensity  scattered  by  the  sphere  in  anv  direction,  assuming 
unpolarized  incident  light,  is  the  average  over  all  possible  orienta- 
tions of  the  incident  electric  vector  or 


I (m,x,8)  - -iL  (i  + i ] 

2kZ  1 ' 1 

in  which  k is  the  propagation  constant,  k = 2^/^,  and  I0  is  the 
incident  flux. 
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The  total  extinction  and  scattering  cross  sections  and  a can 

also  be  expressed  in  terms  of  the  Mle  coefficients 
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The  notation  Re  denotes  the  real  part  of  the  complex  quant  1 tv  In  the 
brackets  , while  and  Qs  are  the  extinction  and  scattering 

efficiencies  defined  as  the  corresponding  cross  sections  normalized 
by  the  cross  sectional  area  of  the  sphere. 


The  true  absorption  efficiency  (cross  section)  is  defined  to  be 
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Thus  for  nonabsorbing  dielectric  spheres,  the  scattering  and  extinc- 
tion coefficients  are  equal.  The  scalar  phase  function  for  a single 
particle  is 


P(Y) 

4ir 
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and  is  obviously  normalized  to  one. 

For  an  element  of  volume  containing  a number  N of  identical 
spherical  particles,  the  volume  scattering  across  section  is  simply 


0 (m,x)  » a N ■ NO 

s 3 4 s 
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vith  dimensions  of  cm  /cm  or  cm  . The  differential  volume 
scattering  coefficient  is 

V 

; (n,x,8)  = N t (m.x,1:)  = f i (m,x,h)  + i (m,x,0)] 

s s 2fr  2 1 

arc  the  scalar  phase  function  for  the  volume  element  is 


PCY)  _ P s ( m , x , ^ ) 

4-  s(n,x) 

= [i7(m,x,-;)  *•  i,<m,x,e)] 

2-x^Qg(m,x) 

ich.  is  identical  to  the  phase  function  for  a single  particle.  The 
phase  function  is  independent  of  the  amount  of  scattering  material  in 
tve  volume  and  depends  only  on  the  particle  size  parameter  x and  the 
refractive  index  in  (assuming  that  the  volume  does  not  contain  so 
mar."  particles  that  multiple  scatte-ir t -refects  appear). 

As  one  might  imagine,  the  scalar  phase  function  plays  a significant 
role  in  single  (and  also  multiple)  scattering  problems.  In  fact,  much 
vrrV.  has  been  done  to  determine  under  ••■•hat  conditions  the  entire 
Vie  series  is  not  required  and/or  can  be  replaced  by  some  analytical 
expression  [5,9,12,13].  We  will  find  evidence  of  this  quest  in  a 
later  subsection  concerned  with  puls«-'  propagation  through  a multiple 
scattering  medium. 

2.i.:  elements  of  Radiative  Transfer  Theorv 

"Vice  the  relevant  single  scattering  properties  are  known,  it  is 
p :s  s '!> 1 e to  •••rite  iown  the  emission  confident  explicit  lv.  Consider 
a cvlindrical  element  of  volume  (f  igu  re  2)  characterized  by  a volume 
-.■•.traction  coefficient  Pext,  a scalar  phase  function,  and  an  albedo 
:r  single  scattering  Assuming  the  volume  has  a base  area  dA  a^a 

--  'eight  dz,  consider  a beam  of  intensity  I contained  in  a solid  angle 
c~'  that  is  incident  on  the  top  surface  of  the  volume  at  an  angle  6’ 

: tee  surface  normal  and  at  an  arimuthal  angle  8'  from  a reference 

plane  through  the  normal.  the  amoun'  of  energy  entering  the  volume  ir. 
it  time  is 


T (8  ',<!>'  1 d^j'  dA  cos  8' 
s-.i  the  path  through  the  volume  fs 
dr  sec  0'  . 


9 


The  energy  lost  from  the  beam  in  unit  time  is 
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or 


6 I(0',d>')dA  cos  B'du'dz  sec  0' 

Kext 


Pext  * ( 6'  . <t>'  )<lui'dAdz  . 


Of  this  energy,  a fraction  cj0  will  be  lost  due  to  scattering  bv  the 
volume.  The  scattered  energv  will  be  distributed  over  directions 
(0,4>)  into  solid  angle  du  according  to  the  phase  function  P (\ ) where 
V is  the  angle  between  directions  ( © * , 4>* ) and  (0,4>).  lae  energv  per- 
unit  time  interval  that  is  scattered  into  the  solid  angle  duj  from  the 
volume  element  dV  = dAdz  will  be 

P(Y) 

(3  I ( 0 ' , <}> 1 ) dojduj ' . 

next  4 rr 


The  emission  coefficient  js  is  defined  as  the  a-runf  of  energv  per 
unit  time  interval  emitted  by  an  element  of  volume  in  a given  direction. 
By  integrating  over  all  possible  incident  angles  ( 0'  we  have 

fw  r 

j (0,4>)  dV  dojdt  = w dV  dwdt  — — / P C>)  1 1©  ’ ,4>')  dw'.  (7) 

® O 4TT 

The  source  function  for  an  elemental  volume  is  defined  to  be 


P r’(G,<t>)  = js  ( 8 ,<±> ) . 
ext 
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Let  us  now  use  these  latter  two  quantities  to  solve  the  equation 
of  transfer  in  more  detail.  Figure  3 depicts  the  problem's  geometry. 
A plane-parallel  homogeneous  and  isotropic  medium  is  illuminated  from 
a direction  (8 0 ,cp 0)  relative  to  the  surface  normal  bv  plane-wave 
irradiance  of  r F watts/cm  . Since  we  are  assuming  the  medium  to 
be  homogeneous  and  isotropic,  the  single  scattering  properties  are 
independent  oc  position  inside  the  medium,  but  the  intensity,  i.e., 
source  function,  depends  on  depth  of  penetration  and  on  direction. 

It  is  convenient  to  use  a normalized  optical  thickness  defined  as 


dt  = p d-, 

ext 

as  the  depth  variable.  Thus,  along  anv  incremental  path  ds  inside  the 
slab  structure  making  an  angle  G to  the  normal , we  have 


ds  = dz  sec  0 = sec  0 . 


11 


This  implies  we  can  write  the  equation  of  transfer  as 

di(t,e,4>) 

cos  e — *=  - i(T,e,  40  + 8(t,e,4>). 

dr 

The  intensity  field  I(t,6,4>),  as  you  recall,  is  composed 
of  two  distinct  parts:  the  incident  intensity,  attenuated  by  the 
extinction  processes  incurred  along  the  path  of  propagation,  and  a 
diffuse  intensity  due  to  scattering  and  true  emission.  Here,  we  again 
chose  to  ignore  true  emission,  and  consider  only  scattering-based 
emission.  Therefore,  let  us  write  the  source  function  as  the  sum  of 
two  terms 


B(r,e,4>)  = Bs(i,6,4>)  + Bo(t,6,4>), 


where  Bs(t,6,4>)  denote  the  source  function  due  to  scattering  into 
direction  (0,4>)  of  the  diffuse  intensity  field  incident  on  the  volume 
element  located  at  optical  depth  t.  Bo(t,0,4>)  is  the  source  function 
due  to  the  scattering  of  the  incident  radiation  by  the  volume  eTenent 
from  the  direction  (6  ,4>  ) into  (0,4>) . Mathematically,  it  follows 
from  equations  (6)  ana  (9)  that 


Bs(-,0,<b' 


uo 

4tt 


2n  tt 


P(v)T(T,e,4>')  sin  e'de'd^' 


o o 


where  I(  r,0,4>)  denotes  the  diffused  intensity. 

The  direction  source  function  Bo(T,6>40  comes  above  as  follows: 
Initially  it  was  assumed  that  the  incident  irradiance  of  ttF  watts/cm 
illuminated  the  planar  structure  at  angle  (80,(J>0).  The  irradiance 
reaching  an  optical  depth  r along  a path  inclined  at  an  angle  0o  to  the 
normal  axis  of  the  structure  is 


ttF  exp  J -t  sec  0q  ' . 


Thus  we  have 


P<V 


' O'  i 

Bo(t,0,4>)  = wo  — ■j — F exp  j - t sec  0^ 


(R) 


j 
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with  y being  the  angle  between  (®o»^o^  an<^  This  implies 


B(T,e,4>) 


<T,e' 


4>')  P(y)  sin  e'  d6'd4> 


t sec  6 

ot 


(<» 


where  P(y)  = P(0 ,$ ,9 ' A ' ) . One  need  only  apply  the  appropriate 
boundary  condition  in  order  to  complete  the  solution  to  the  problem. 

2 . 3 Mutual  Coherence  Function 

In  the  above  analysis,  we  regarded  radiation  transfer  purely  as  a 
collection  of  photons  diffusing  through  a scattering  medium.  This  is 
sometimes  a dangerous  oversimplification  of  the  problem  as  one  has 
neglected  to  investigate  the  evolution  of  the  radiation  field's 
coherence.  One  finds  in  practice  that  signal  coherence  can  affect 
the  performance  of  communication  system  significantly,  e.g.,  a loss 
of  spatial  and  temporal  coherence  in  a heterodyne  system  can  cause 
channel  fading.  Thus  one  recognizes  propagation  through  single-  and 
multiple-scattering  media  to  be  a stochastic  process. 

Stochastic  processes  are  generally  analyzed  in  terms  of  their 
statistical  moments  [14,15].  Probablv  the  most  important  statistical 
moment  is  the  second  order  moment;  this  is  because  of  the  relative  ease 
in  which  it  can  be  handled  both  mathematically  and  experimentally.  One 
finds  that  many  of  the  theoretical  studies  concerned  with  scattering 
media  utilized  this  technique.  For  example,  Heggestad  [19]  derived  a 
generalized  scattering  function  from  this  moment  for  his  paper  on 
optical  communications  through  multiple-scattering  media.  Ishimara 
[14],  on  the  other  hand,  used  the  two  frequency  mutual  coherence 
function  (second  order  statistical  moment  of  the  field)  in  his 
analysis.  Livingston  [17]  and  Lutomirski  [18]  were  both  successful 
in  developing  a generalized  mutual  coherence  analysis  of  multiple- 
scattering effects  on  optical  propagation.  Since  all  of  the  papers 
are  based  on  the  same  relative  formulation,  we  will  describe  only 
Livingston's  work  and  leave  it  to  the  reader  to  investigate  the  other 
variations  cited. 

2.3.1  Derivation  of  the  Radiation  Transport  Equation 

Consider  a semi-infinite  slab  of  scattering  medium  with  faces 
perpendicular  to  the  axis  of  a coordinate  system.  Let  the  slab  have 
finite  thickness  in  the  7.  direction.  The  radiant  intensity  S(r,n)  is 
conventionally  defined  as  the  flux  of  energy  through  a unit  area  in 
unit  time  dt  moving  in  the  mean  direction  denoted  bv  the  direction 
vector  n.  That  is 

dE  joules  ■ dA  S(r,n)  dndt  (10) 
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or 


S(r,n)  “ dE/dAdndt  watt s/m^  • ater 


(ID 


where  dA  is  an  element  of  area  and 

dri  * sin  n d(di  (12) 

is  an  element  of  solid  angle.  Equation  11  defines  a pencil  of  radiation 
Introduced  by  Chandrasekhar  [8].  As  an  elementary  area  of  a wavefront 
dA  moves  with  the  wave  at  velocity  c,  it  sweeps  out  a volume  as  shown 
in  fig  4.  During  this  passage  radiant  energy  may  be  gained  or  lost, 
depending  on  whether  energy  is  scattering  into  the  volume  moving 
parallel  to  n or  lost  to  it  by  extinction.  The  radiative  transfer 
equation  then  represents  a balance  of  gains  and  losses  during  this 
process.  For  the  moment,  we  will  suppress  the  dependence  of  the  cross 
sections  on  particle  size  and  refractive  index.  For  polydisperse 
systems  average  radiation  transfer  is  obtained  by  summing  over  all 
particle  sizes.  Let  n be  the  number  of  scattering  centers  per  cubic 
meter  and  aext,  asca  and  aa{,s  be  the  total  extinction  cross  section, 
total  scattering  cross  section,  and  total  absorption  cross  section, 
respectively,  defined  in  the  previous  section.  Assuming  no  shadow 
effects,  the  fraction  of  the  elementary  area  covered  by  scatterers  is 

A. 

fraction  ■ iext  n(h*Z)c  dt.  (13) 


Therefore  the  loss  to  the  radiant  energy  due  to  extinction  is  given  by 
multiplying  Eq.  13  by  S dAdndt  to  give 

d2ElogH  = foext  n(ri*Z)c  dt](S  dAdndt).  (14) 

Similarly,  radiation  scattered  into  the  volume  moving  in  the  direction 
n is  given  by 

d~E  = o n(n'Z)c  dt  ; — f P(h ,n ' ) S(r ,n ' ) dri 'dAdt . (15) 

gain  sc  l 

Here,  P is  the  probability  of  scattering  radiation  from  direction  n'  to 
directior  n.  The  net  change,  gain  minus  loss,  is  just  °qual  to  the 
substantial  derivative  of  E,  namely,  (dF./dt)dt.  For  steadv-state 
conditions 


dt 


• VEdt  = cn’VS  dAdh(n*Z)dt, 


(16) 


where  the  defining  equation  for  S has  been  used  and  also  the  fact  that 
the  wave  moves  with  the  speed  of  light  c.  Thus  the  equation  satisfied 
by  S,  namely 


no 


?S(r,n) 


Jext 


nS(r ,n) 

'Xj 


4tt 


J P(n,n’ 


)S(r,*,)dn’, 


(17) 
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is  the  standard  radiation  transport  expression.  In  the  preceding,  a 
slightly  more  general  expression  follows  if  we  allov?  S to  denend  on 
the  light  frequency.  u>,  as  would  be  required  if  Doppler-scattering  is 
of  interest: 


n • VS(£,n,iu)  + naext  S(r,n,w)  - 


noSc 

4rr 


w'iSCr.n.w^dn' . 


(1R) 


Apart  from  differences  in  notation  and  a more  complex  kernel,  this 
equation  is  identical  with  Fante's  [19].  Either  Eq.  17  or  Eq . 18  along 
with  an  expression  for  P serve  to  define  the  evolution  of  the  radiant 
intensity  S through  the  cloud.  Since  P can  be  expressed  in  a series  of 
Legendre  polynomials,  Eq.  18  can  be  numerically  treated  as  an  infinite 
set  of  coupled  differential  equations  or,  what  is  equivalent,  a set  of 
homogeneous  algebraid  equations  to  be  solved  by  matrix  inversion. 

Before  considering  the  relationship  of  S to  quantities  of  interest, 
let  us  rewrite  Eq.  17  as  follows: 


3 

3T 


S(p,T,n) 


+ 


n L 


naext 

^o  r 

J 


•VI  S(p,T,njL)  + S(p,T,n 
dniJP(ni.,ni.  )S(p,x  ,n'x)  , 


1) 


where  the  optical  depth  t is  given  by 


T = naext  Z 


(19) 


and  the  albedo  for  scattering 
%m  °s/°ext  S 1 


has  been  introduced.  In  the  notation  of  Eq.  19,  nj_ is  the  perpendicular 
direction  to  Z.  Further  progress  in  analyzing  Eq.  19  can  be  made  if  S 
depends  on  t and  the  central  angle  alone,  that  is,  if  S is  homogeneous  in 
planes  transverse  to  Z,  or  £ depends  only  on  the  difference  of  the  vectors 
nl  and  n i.‘  The  former  case  is  fully  considered  in  Chandrasekhar's 
work,  and  the  latter,  termed  .small  angle  scattering  approximation, 
is  considered  by  Fante  in  Ref ,17 . Let  us  define' 


F(k,£,T)  -(2*) 


illl 


S(^,nL,i)e 


n L 


drdn 


(20a) 


16 


and 


Q(p)  - (2ir)' 


d?lP(5lL)e 


ilW 


(20h) 


where  it  is  now  assumed  that 


P(ni  ,n'j_)  = P(nj_  - n'L  ) . 


(20c) 


The  equation  satisfied  by  F is  a combination  of  Fqs.  10,  20a,  and  20b,  and 
the  result  is 

~ (k , p , t)  + F(k  ,p,  t)  - F(k  , p , t)  - TTUQ(p)F(k,p,T)  = 0.  (21) 

d t % V nCT  'V,  % U rU 


The  solution  of  the  first-order  equation  in  the  optical  depth  is 
readily  effected: 

F(k,p,  x)  - F(k,p,0)e-  T 

The  inverse  transform  of  this  equation  yields  the  solution  for  S: 


S(r  ,n±,r  ) = e 


7 //. f-  ’)  <«, 


Here,  F'(k,p',0)  is  the  transform  of  the  initial  radiance  distribution, 
and  it  remains  only  to  specify  0(p)  to  obtain  an  explicit  solution. 
Oiven  an  initial  Held  n(r),  then,  the  initial  transformed  irradiance 
distribution  can  be  written 


:,£,0)  = (P7!)  j e'x  * u(r+p/2ko)u*(r-p/ 


2k  ) dr, 

O 'u 


where  kQ  is  2rr />, . A similar  expression  holds  for  anv  value  of  t and, 
when  averaged  over  a suitable  ensemble,  is  related  to  the  mutual 
coherence.  Thus 


H 

II 


The  last  expression  follows  after  the  substitution  of  the  solution, 

Eq.  22,  into  the  second  line.  It  is  a mutual  coherence  function,  defined 
in  Eq.  25,  that  is  central  to  both  a description  of  the  propagation 
process  and  to  the  carrier-to-noise  calculation  of  the  optical  heterodyne 
communication  system.  In  Eq.  25,  R is  the  center-of-mass  vector  in  the 
receiver  plane  (perpendicular  to  tfre  axis  of  propagation)  and  p is  the 
difference  vector  in  the  same  plane: 

lm  2 + 


and 


- r0 


where  r3  and  r^  are  two-dimensional  vectors  in  the  receiver  plane.  Let 
the  lauhched  beam  have  the  form 


„(,)  - + (‘X0r2/2Eo)  (26) 

where  A is  the  effective  aperture  area  and  RoIs  the  radius  of  curvature 
of  the  faunched  wavefront.  (A  collimated  beam  is  g’ven  hv  F =*  -■*>.) 
Combining  Eqs . 22  and  2D  leads  to 


kap+ 
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Ao|uo!  -(Pik2+62P2+63k-p) 
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where 
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Substitution  of  Eq.  27  into  the  expression  for  T leads  to 

a I |2  ** 

Aq  I Uo  I 

r (r,p,t)  r 

^ 4ir 


) the  expression  tor  i leans  to 
^-t-gjk^p2  j~  81k2-63k*p-ik-R^>u  9(kcp)' 
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The  integral  over  k can  be  carried  out  to  give 


2 2 


r(R,p,x) 


A lu  I 
o o 

4*0, 


*.p/2.  , . (30, 


For  a collimated  beam  centered  at  R = 0 

2/  / / 1 /_ * 


A ' 12 

A i u 


o o 


r(0*E’T)  = 


2 2 

, w z 

A + — ~ — 
o . 2 


T-[p2/4(Ao/TT+TrZz/koA0)]-^-uQ(koP)T  . (31a) 
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If  we  further  assume  that  the  receiver  is  in  the  far  field  of  the 
transmitter  aperture,  we  can  write  out  a simple  expression  for  the 
normalized  mutual  coherence  function: 


y(O,0,t)  - + •».»' 


H2a) 


2.4  Monte  Carlo  Computer  Simulation 

Monte  Carlo  computer  simulation  can  be  applied  to  any  problem  if  one 
knows  the  probability  for  each  step  in  a sequence  of  events  and  desires 
the  probability  of  the  total  of  all  possible  events.  Many  researchers, 
most  notably  Bucher  [12],  and  Plass  and  Kattawar  [18],  have  contributed 
significantly  to  the  area  of  laser  irradiance  through  single-  and  multiple- 
scattering media  using  this  technique.  The  advantages  of  this  method  are 
the  following  [20]: 

(1)  The  calculation  may  be  performed  for  any  single  scattering  function 

regardless  of  the  degree  of  anisotropy.  j 

(2)  Any  relevant  parameter  such  as  the  single  scattering  albedo,  number 
density  of  molecular  and  aerosol  scattering  centers,  and  the  various  cross 
sections  may  be  varied  with  optical  depth  in  the  atmosphere  in  any  desired 
fashion. 

(3)  The  radiation  reflected  by  the  boundary  interface  may  follow  any 

desired  distribution.  j 


(4)  Any  reasonable  number  of  polar  and  azimuthal  angles  may  be 
selected . 

(5)  The  average  mean  optical  path  of  both  the  reflected  and  trans- 
mitted photon  may  be  easily  calculated. 
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The  disadvantage  of  the  Monte  Carlo  method  is  twofold:  First,  the 

standard  deviation  of  the  results  is  roughly  inversely  proportional 
to  the  square  root  of  the  computing  times*  Therefore  this  technique 
may  not  be  practical  if  a high  degree  of  accuracy  is  required.  Second, 
this  technique  does  not  lend  insight  into  the  physics  of  the  process  nor 
does  it  immediately  provide  the  functional  dependencies  required  for 
parameter  variation. 

The  purpose  of  this  section  is  to  review  the  Monte  Carlo  method 
in  more  detail  utilizing  a representative  paper,  reference  [12]  by  E . A. 
Bucher  to  be  exact. 

2.4.1  The  Method 

A Monte  Carlo  Computer  simulation  is  essentially  a calculation 
based  on  some  mathematical  model  of  a physical  situation.  We  will 
now  describe  the  model  and  technique  Bucher  employed  in  his  analysis. 

The  configuration  chosen  for  the  multiple  scattering  medium 
was  a plane-parallel  slab  of  infinite  width  and  physical  thickness  T. 

The  medium  Itself  was  assumed  to  be  homogeneous  and  isotropic  in  nature. 
The  above  structure  was  chosen  for  its  mathematical  simplicity. 

Figure  5 depicts  a typical  propagation  path  simulation.  A normally 
incident  ray  travels  a distance  d^  until  it  is  scattered  the  first  time. 

On  scattering  a new  direction  (0^,  4^)  is  selected  relative  to  the 
incident  ray.  The  ray  then  travels  a distance  d2  to  its  second 
scattering.  At  each  successive  scattering,  a new  direction  (e^,  <J>^) 
relative  to  the  incident  normal  is  chosen.  Then,  the  ray  travels  a 
distance  d1+^  to  the  next  scattering  event.  This  process  continues 
until  the  ray  leaves  the  scattering  medium  bv  either  of  its  two  planar 
boundaries.  Ravs  leaving  through  the  top  (incident)  boundary  are 
indicative  of  back-scattered  irsadiance.  Rays  exiting  the  bottom  boundary 
indicate  the  transmitted  light.  For  each  of  these  latter  ravs , the 
computer  program  records  its  total  travel  distance,  exit  angle,  and 
radius  of  the  exit  point  relative  to  the  straight-through  axis.  This 
information  was  used  later  to  calculate  several  statistical  measures 
after  a sufficient  number  of  propagation  paths  had  been  traced. 

The  distance  dj  that  a ray  travels  in  going  from  the  (.]-!) 
scatterer  to  the  I**1  scatterer  was  set  up  by  Bucher  to  be  randomly 
selected  from  a probability  density 


P(dj)  - (-^)  exp  {-  d./D}  ; 

J D J 

That  is,  the  random  variable  d was  chosen  to  assume  values  between 
d and  d+Ad  with  total  probability  (1/D)  exp  {-  d/D)  Ad.  The  quantity  D 
is  the  scaling  parameter  of  the  exponential  and  is  also  the  average  value 
of  the  random  variable  d. 
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The  scattering  angle  p.  and  are  selected  Independently  for  each 
scattering  event.  The  angles  9 and  $>  represent  the  polar  and  rotation 
angles,  respectively,  in  a spherical  coordinate  system  whose  origin 
is  in  the  direction  of  the  Incident  ray.  In  this  configuration,  9 is  the 
angle  between  the  incident  and  scattered  rays  in  the  plane  of  scattering: 
and  p is  the  rotation  of  the  scattering  plane  about  the  incident  ray. 

At  each  scattering,  <p  was  randomly  selected  from  all  possible  values  of 
5 between  0 and  Zv. . This  uniform  distribution  for  the  rotation  angle 
can  be  justified  either  on  the  basis  of  spherical  scatterers,  or  on 
tiie  basis  of  arbitrarily  shaped  scatterers,  with  no  preferred  spatial 
orientation . 

Selection  of  the  0-angle  distribution,  i.e.,  scalar  phase  function, 
involves  the  particular  scattering  mechanism  of  the  medium  and  raises  both 
physical  and  computational  questions.  As  noted  before,  we  will  return 
to  this  question  in  the  section  on  pulse  propagation  through  mul.tinle 
scattering  media. 

2.4.2  The  9-angle  Distribution  Function 


During  our  discussion  of  Monte  Carlo  computer  simulation,  we 
stated  that  the  choice  of  the  9 -angle  distribution  function  involved 
the  particular  scattering  mechanism  of  the  medium  modeled  and  thus 
raised  both  physical  and  computational  questions.  Bucher  responded 
to  this  situation  by  using  two  approximate  scattering  functions  in 
conjunction  with  the  full  Mie  series  scattering  function  in  his 
analysis  and  then  comparing  their  results.  The  two  approximating 
functions  chosen  were  the  megaphone  and  the  Henyey-Creenstein 
functions.  As  it  turned  out,  Bucher  found  good  agreement  between  all 
three  scattering  functions  for  diffusion  thicknesses  greater  than  3. 
The  diffusion  thickness  is  a weighted  optical  thickness  and  is 
lefined  to  be 

■ , * (1  - < cos  8>  ) T/D 


where 

T 3 physical  thickness  of  the  medium 

D £ usual  mean-free  path  between  scatterings 

<cos9>  - average  cosine  of  the  scattering  angle  relative  to 
the  incident  d^-ection 

.he  wavelength  used  In  the  analysis  was  .7p  and  the  particles  were 
assumed  to  have  a refractive  index  of  1.34  and  a radius  distribution 
'-etween  3.8y  and  4 . S u . 
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2.5  Theoretical  Results 


2.5.1  Multipath  Propagation  Parameters  [4,12,16,20-23] 

If  a collimated  optical  pulse  Is  Incident  on  a multiple-scattering 
region,  the  beam  will  experience  beam  spreading,  dispersion  in  angle 
of  arrival,  attenuation,  degradation  of  spatial  coherence,  and  dispersion 
in  time  and  frequency  of  the  signal  modulating  it  after  transversing  the 
region.  These  effects  are  related  to  the  following  multipath  propagation 
parameters . 

a.  angular  spreading 

b.  spatial  spreading 

c.  multipath  time  spreading 

d.  total  transmission 

e.  Doppler  spreading 


Figure  6 illustrates  a few  of  these  parameters  graphically.  We  have 
omitted  polarization  spreading  as  Bucher  chose  not  to  include  it  in 
the  simulation. 

"Angular  spreading"  constitutes  an  effective  decollimation  of  the 
incident  radiation.  It  contributes  to  beam  spread,  dispersion  in 
angle-of-arrival , and  loss  of  spatial  coherence. 

"Spatial  spread,"  on  the  other  hand,  indicates  the  dimensional 
increase  of  the  beam's  finite  cross  section.  It  is  related  to  the  above 
parameter  In  that  angular  spread  inside  the  medium  produces  a spatial 
spreading  as  the  pulse  propagates  on.  Spatial  coherence  and  beam 
spread  are  related  to  this  parameter. 

Different  distances  along  the  various  possible  propagation  paths 
imply  different  transit  times  for  a photon.  Thus,  a short  optical 
pulse  will  incur  pulse  broadening  after  transversing  a multiple 
scattering  region.  Communication  theorists  often  refer  to  this  effect 
as  "multipath  time  spreading"  and  denote  its  counterpart  parameter  by 
the  letter  L and  the  same  name  [21,25].  One  finds  that  this  quantity 
is  very  important  in  communication  applications  in  that  1/L  represents 
the  maximum  pulse  frequency  which  can  be  resolved. 

"Total  transmission"  describes  the  amount  of  irradiance  left  in  the 
pulse  after  transversing  the  medium.  It  is  directly  related  to  the 
beam  attenuation. 

If  the  scatters  in  the  medium  are  moving,  one  finds  that  a 
change  in  the  wavelength  (frequency)  of  the  light  will  occur  after 
each  scattering  event.  This  wavelength  change  is  called  a "Doppler 
shift."  If  these  scatters  are  also  in  relative  motion  to  one 
another,  the  total  Doppler  shift  incurred  along  the  various  scatter 
paths  will  be  different.  Hence,  the  irradiance  transmitted  through 
the  scattering  region  from  a single  frequency  source  will  contain  a 
range  of  "Doppler  spread"  frequencies.  We  will  ignore  this  parameter 
in  the  present  discussion  since  current  optical  receivers  fail  to 
possess  both  the  selectivity  and  the  sensitivity  necessary  to  measure 
this  effect. 
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For  simulations  involving  clouds  with  t,  2,  Bucher  found  the 
distribution  of  exit  angles  relative  to  the  cloud's  normal  to  be 
fairlv  uniform  in  brightness.  That  is,  he  observed  that  the  exit 
angles  for  ravs  emerging  from  a given  area  on  the  bottom  boundary 
tended  to  be  concentrated  within  a 45°  cone  about  the  normal. 

However,  when  a projected  area  correction  factor  of  1/cos  <1^,  4- 
being  the  exit  angle,  was  applied  to  the  data,  the  resulting 
distribution  did  not  contradict  a hypothesis  of  uniform  brightness 
in  the  bottom  boundary  when  viewed  from  all  angles. 

Vith  thin  simulated  clouds,  Bucher  observed  a sharply  peak 
distribution  of  exit  angles  centered  around  the  straight-through 
axis.  In  addition,  the  peak  was  found  to  contain  much  more  trans- 
mitted radiation  than  just  the  contribution  of  the  direct  extinction 
beam.  Ke  attributed  this  additional  irradiance  to  a strong  forward 
scattering  behavior.  In  particular,  he  found,  using  the  Mie  series 
scattering  function,  that  the  peak  was  most  pronounced  within  5° 
of  the  normal  and  that  the  contribution  was  generated  within  0.08T 
of  the  central  point  where  the  unscattered  ray  would  exit.  Figure  7 
shows  a plot  of  the  fraction  of  total  incident  rays  exiting  the 
cloud  generated  during  the  simulation  as  a function  of  optical 
thickness  in  relation  n-  the  unseat terod  component.  Bucher  postu- 
lated that  this  multiple  forward  scatter  component  may  be  quite  useful 
in  some  optical  comnunicat ion  situations,  o.g.,  over-t he-hor izon 
eonr.unicat  ions . 

2.5.3  Spatial  Spreading 

Scat  ial  spreading  was  measured  bv  Bn  .her  in  terms  of  two  criteria. 
Both  of  these  criteria  were  based  on  the  radius  r defined  previously 
as  the  distance  from  the  exit  point  to  the  straight-through  axis  on 
the  '-ottor  boundary  (see  fig  51.  ’’v"  • all'’,  the  r distribution,  had 

lone  tails  with  some  rays  having  a rad lus  several  times  the  average. 
Thus  Bucher  chose  to  characterize  spatial  spreading  bv  the  average 
-i1;  - r and  the  exit  radius  r inside  which  half  the  transmitted 
rav-  exited  the  medium.  Figure  8 shows  <r>/D(j  as  a function  of 
for  the  various  scattering  functions  used.  Figure  **  shows  rc/R^  as 
: function  of  j . The  parameter  9^  is  given  bv  the  equation 

= n/(l-<cosu>) 


where  both  9 and  -cos  ' have  boor,  defined  previously.  Figure  10 
depicts  <cos  as  a function  of  the  normalized  Mie  scattering  parameter. 
For  the  scattering  medium  modeled,  tD/A  ranged  from  17.05  to  .-1.54,  therebt 


value  of  .81  to: 


b *■  • *-  '*6 

apparent  from  figures  7 and 
functions  occur  for  tj  N3. 
mate  these  curves  using  the  equations 


cos  This  implies  t ^ - 5.26  T.  It  is 

8 that  common  behavior  for  all  three  scattering 
Furthermore,  Bucher  found  that  he  could  approxi- 
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<r>  = 0.92  0 (id)' 
r c = °‘78  nd(Td)( 


(33a) 

(33b) 


respectively,  for  rd>3.  Recalling  the  relationship  between  T and  0d , we 
see  that  the  above  expressions  can  be  rewritten  as 


<r  = 0.92  T (Td) 


r c = 0.78  T (rd) 


-0.08 

-0.07 


(34a) 

(34b) 


respectively.  With  this  normalization,  we  see  that  the  spatial 
spreading  saturates  and  even  decreases  as  a cloud  of  fixed  thickness 
becomes  more  opticallv  dense.  This  is  reflected  in  figures  11  and  12 
for  <r>  and  r , respectively. 

2.5.4  Multipath  Time  Spreading 


In  character izing  multipath  time  spreading,  Bucher  found  it 
necessary  to  use  three  statistical  measures:  the  mean  pulse  width  <f, 
its  standard  deviation  a , and  the  width  of  the  most  intense  central 
portion  of  the  time  distribution.  This  latter  parameter,  which  is 
the  multipath  time  spread  L,  is  defined  as  the  shortest  range  of  multi- 
path  time  values  which  encompass  63%  of  the  total  transmitted  rays. 

pigures  13-15  show  ct  , <t>,  and  L,  respect ivelv , as  a function  of 
the  diffusion  depth  ;d  for  a variety  of  optical  thicknesses  ranging 
from  5 to  120.  It  is  apparent  that  all  three  scattering  functions 
give  approximatelv  the  same  measures  for  values  of  Td>2.  The  normal- 
ization factor  tj  used  in  these  figures  is  given  bv  the  equation 

t = O/c 
d d 

where  c is  the  speed  of  light  in  air.  Bucher  then  went  on  to  show 
that  the  three  statistical  measures  could  be  presented  as 


<t  > = 

0.62 

fcd 

(35a) 

"t 

0.64 

ld 

(Xd)1-81 

(35b) 

L 

0.74 

t-d 

(Td)1’88 

(35c) 

for  td >3 . Tt  can  be  easily  shown  that  the  above  expressions  can  be 
rewritten  as 
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It  is  easily  seen  from  the  above  that  the  time  spread  for  sufficiently 
thick  clouds  grows  in  a way  that  is  essentially  independent  of  both 
the  details  of  the  scattering  function  and  of  the  width-measuring 
criteria . 

2.5.5  Total  Transmission 

It  can  be  shown  [22,23]  that  the  fraction  of  incident  rays  trans- 
mitted through  a "thick"  cloud  of  nonabsorbing  scatterers  is  given  by 

It  - A /(id  + B),  (37) 

where  A is  a function  of  the  incidence  angle,  and  B is  essentially 
a constant.  For  normally  incident  radiation,  this  asymptotic 
expression  becomes 

It  = 1.69/(xd  + 1.42) . (38) 

Figure  16  shows  a plot  of  simulation  Ifc  and  the  It  predicted  from  the 
above  equation  as  a function  of  t<j.  As  one  can  see,  there  is  very  good 
agreement  between  the  two  formulations.  Thus  we  have  an  analytical 
expression  for  It  in  terms  of  the  weighted  optical  thickness. 

3.0  EXPERIMENTAL  TEST  PLAN 

3.1  Goals  of  the  Program 

The  basic  measurement  proposed  here  is  to  determine  the  time  history 
of  a short  laser  pulse  that  has  propagated  through  a path  of  fog  particles 
characterized  by  drop-size  distributions  and  spatial  extent.  The  power 
loss,  the  pulse  spreading  and  the  pulse  delay  will  be  measured  as  a 
function  of  receiver  FOV  for  various  types  of  marine  coastal  fog.  Direct 
beam  extinction  coefficients  will  also  be  measured  and  compared  with  that 
derived  from  Mie  theory  using  the  measured  drop-size  distributions  and 
readings  from  a visiometer.  The  laser  beam  for  the  experiments  will  be 
supplied  by  a frequency  doubled  (.53  um)  neodymium  YAG  laser  with  an 
automated  digital  data  acquisition  system.  The  fog  particle  size  distribu- 
tions will  be  measured  using  two  Knollenberg  drop-size  spectrometers 
situated  selectivelv  along  the  propagation  path.  For  fog  events  which 
occur  during  daylight  hours,  their  time  history  and  spatial  extent  will  be 
determined  from  time-lapse  photographs  taken  at  two  separated  locations. 
During  night  hours  other  atmospheric  sensors  (discussed  in  a following 
subsection)  will  be  utilized.  The  laser,  variable  FOV  receiver,  drop- 
size  spectrometers  and  other  meteorological  sensors  to  be  used  in  the 
measurements  are  described  in  detail  in  different  sections  which  follow. 
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Figure  16.  Total  transmission  as  a function  of  r^.  The  theoretical 
curve  is  from  van  de  Hulst's  work.  The  909f  confidence  limits  on 
all  data  points  are  ±59 < or  better.  [12]. 


3.2  Rationale  for  Selection  of  Site 


Ideally,  the  measurements  should  be  conducted  from  an  aircraft  to  a 
receiver  on  the  ground.  However,  the  cost  and  equipment  available  presently 
prohibit  this  approach.  An  Equally  acceptable  configuration  would  be 
through  a vertical  section  of  clouds  down  a mountain  slope  or  between  two 
mountain  peaks  where  the  clouds  are  typically  of  marine  origin  (as  exists 
on  the  island  of  Hawaii) . However,  for  the  preliminary  measurements  a 
site  located  on  Point  Loma,  San  Diego,  California,  has  been  chosen. 

This  site  is  a compromise  based  on  economics  and  time.  The  selected 
propagation  path  is  shown  in  the  aerial  photograph  of  Point  Loma,  figure 
17.  The  two  terminals  of  the  propagation  link  are  separated  by  2.4  km  and 
located  at  NELC's  Battery  Ashbum  facility  and  the  remote  sensing  facility. 
Building  T323.  The  elevations  above  MSL  are  377  feet  and  105  feet,  respec- 
tively. 

The  advantages  of  the  Point  Loma  location  are:  (1)  easily  accessible 
and  located  on  government  property;  (2)  the  propagation  path  is  in  proximity 
to  the  remote  sensing  facility  for  assessing  the  propagation  medium; 

(3)  considerable  marine  fog  (visibility  less  than  1 mile)  occurs  along  the 
coast  during  the  fall  and  winter  months  (see  figure  18) ; and  (4)  other 
conditions  of  reduced  visibility  (less  than  5 miles)  also  frequently  occur 
in  the  winter  and  summer  months  (see  figures  19  and  20) . The  Fleet 
Weather  Service  located  at  nearby  North  Island  Naval  Air  Station  is  also 
available  for  fog  forecast  such  that  personnel  can  be  alerted  and  equipments 
1 made  operational  during  periods  of  high  probability  for  fog  occurrence. 

Generally,  two  types  of  fog  prevail  along  the  Point  Loma  coast,  namely 
stratus-cloud  and  Santa  Ana  fog.  The  mechanisms  producing  each  fog  type 
appear  to  be  distinctly  different.  Stratus-cloud  related  fog  is  a mid- 
calendar year  phenomenon  and  is  produced  when  the  base  of  the  stratus 
clouds  descends  to  the  ground,  usually  at  night.  Evidence  indicates  that 
a sequence  of  events  initiated  by  radiation  cooling  at  the  cloud  top 
causes  the  cloud  base  to  descend  slowly.  The  visibility  in  this  type 
of  fog  seldom  reduces  below  1/4  mile,  which  is  typical  when  marine 
aerosols  are  involved. 

Fog  related  to  wintertime  Santa  Ana  conditions  appears  to  have  two 
formation  mechanisms.  However,  both  types  of  fog  appear  to  be  transported 
into  the  San  Diego  region  by  the  movement  of  mesoscale  troughs  at  different 
orientation.  One  type  of  Santa  Ana  related  fog  usually  moves  into  San 
Diego  as  a diffuse  or  sharp  fog  bank.  The  other,  and  less  common,  Santa  Ana 
fog  moves  as  a low  thin  deck  (about  30  to  50  m thick)  northward  from  the 
coast  of  Baja  California  into  the  San  Diego  region.  Although  the  fog 
formation  mechanism  and  mesoscale  circulation  associated  with  Santa  Ana 
conditions  are  not  generally  known,  there  appears  to  be  relation  between 
tie  vet  bulb  temperature  of  the  continental  air  which  moves  out  over  the 
ocean  and  the  sea  surface  temperature  which  results  from  ocean  current 
up-welling  along  the  coast.  Because  of  the  smaller-sized  continental 
[ aerosols  involved,  this  type  of  fog  is  much  more  dense  than  the  stratus- 

] cloud  type  and  provides  the  best  opportunity  for  observation  of  the  pulse 

.■5j  stretching  effects.  Also,  since  the  Santa  Ana  fogs  approach  the  propagation 

i path  from  the  west  and  south,  various  conditions  exist  which  allow  the  laser 

| source  and  receiver  to  be  both  or  separately  immersed  in  the  fog. 
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A recognized  disadvantage  of  utilizing  the  Santa  Ana  fogs  for  the 
propagation  experiments  is  that  thev  are  contaminated  hv  continental  aerosols 
and  may  not  be  typical  of  clouds  which  occur  over  remote  areas.  Also,  the 
vertical  gradients  in  the  aerosol  distributions  and  temperature  moisture 
profiles  experienced  in  the  low-level  propagation  path  will  not  be  the  same 
as  for  elevated  oceanic  clouds.  Needed  is  the  ability  to  extrapolate  the 
data  acquired  in  the  coastal  fogs  to  the  open  ocean  environments.  Given  the 
drop-size  distributions  and  spatial  variations,  computer  codes  need  to  be 
developed  to  predict  the  multiple  scattering  propagation  effects  for  various 
types  of  oceanic  fogs  and  clouds.  These  computer  programs  involve  compli- 
cated Monte  Carlo  techniques  which  are  expensive  to  develop  and  process  by 
digital  computer  and  are  beyond  the  scope  of  the  measurement  program  proposed 
here.  However,  the  multiple  scattering  codes  require  as  inputs  the  single 
particle  scattering  Mie  coefficients.  Computer  codes  for  these  Mie  coeffi- 
cients and  their  integration  with  drop-size  distributions  to  obtain  extinc- 
tion coefficients  presently  exist  at  NELC.  These  programs  were  developed 
hv  NELC  and  MEGATEK  Corp . in  support  of  the  NELC  OCCULT  (Optical 
Covert  Communications  Using  Laser  Transceivers)  system.  As  mentioned 
earlier,  the  singly  scattered  direct  beam  extinction  coefficients  will  be 
measured  using  a narrow  FOV  receiver.  These  measured  values  will  then  be 
compared  to  those  calculated  with  the  computer  codes  using  measured  drop- 
size  distributions. 

J . 3 Proposed  Measurements 

As  discussed  in  section  3.1,  the  goals  of  this  experimental  test  plan  are: 

a.  The  measurement  of  the  time  history  of  the  radiant  power  of  a 

light  pulse  after  propagating  through  various  clouds  or  marine  aerosols;  and 

b.  The  determination  and  the  quantification  of  the  various  loss 
mechanisms  inherent  to  the  aerosols. 

The  measurements  and  the  presentation  of  the  data  will  be  in  a form 
compatible  with  the  needs  of  the  optical  communication  systems  engineer. 

That  is,  the  subsequent  results  and  conclusions  derived  from  this  program 
will  enable  the  optical  system  designer  to  predict  communication  link 
performance  and  thereby  allow  system  design  optimization.  Let  us  now 
discuss  the  proposed  measurements  in  more  detail. 

Consider  a laser  source  emitting  Pt  watts  with  beam  divergence 
The  power  received  P^  at  a distance  R with  a receiver  possessing  an 
antenna  of  area  A is  approximately 


P 


R 


(39) 


neglecting  atmospheric  losses.  r f we  now  intersperse  3 turbid  medium, 
absorption  losses  are  negligible  compared  to  the  losses  induced  by 
scattering,  i.e.,  uj  - 1.  The  power  then  received  at  the  receiver, 
assuming  very  thin  clouds,  is  given  by 
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where  kg  t is  the  single  scattering  coefficient  previously  described  in 
section  2.2.1.  With  nominal  values  of  kscat  occurring  in  clouds,  the 
practical  possibility  of  attaining  optical  communication  through  clouds 
using  the  direct  beam  would  be  ruled  out.  For  example,  fair  weather 
cumulus  clouds  possess  a scattering  coefficient  on  the  order  of  5xl0-2  m-1 
If  the  cloud  has  a physical  thickness  of  1 km,  the  path  loss  due  to  single 
particle  scattering  would  be  100  dB . It  is  implied  in  the  direct  beam 
formalism  that  any  photon  scattered  is  completely  lost  and  will  not  be 
collected  by  the  receiver.  This  large  value  of  path  loss,  100  dB,  in 
addition  to  other  incurred  losses  would  severely  inhibit  the  potential  of 
communication  through  clouds.  However,  since  the  forward  scattering 
process  does  not  absorb  photons,  but  merely  redistributes  them  both 
temporally  and  spatially,  collected  scattered  irradiance  could  possibly 
provide  a means  to  enhance  open  channel  communications  availability. 


The  equation  for  received  scattered  power  given  above  (Eq.  40)  assumes 
that  the  scattering  by  atmospheric  clouds  is  isotropic  in  nature,  i.e.,  the 
scalar  phase  function  p(V>’t')  is  defined  to  be  unity.  However,  this 
is  not  the  case  in  real  life.  Calculations  using  Mie  theory  show  that  the 
scattering  by  atmospheric  aerosols  is  not  isotropic,  but  rather  highly 
biased  in  the  forward  direction.  Figure  21  illustrated  a large  scattering 
probability  for  angles  less  then  2 degrees,  and  then  shows  an  order  of 
magnitude  reduction  in  probability  to  a value  which  remains  constant 
for  about  30  degrees.  Beyond  this  angle,  the  scattering  probability 
decreases  by  one  or  two  orders  of  magnitude. 

Figure  22  is  a typical  schematic  representation  of  the  laser  source  as 
it  might  be  viewed  by  the  receiver  after  being  scattered  by  cloud  aerosols. 
The  central  core  (shown  here  as  a dot)  Is  the  unscattered  beam,  which 
possesses  microradian-type  angular  spread.  Surrounding  this  is  an  aurora 
of  width  9^  (on  the  order  of  a few  degrees)  which  consists  of  photons 
which  have  been  scattered  only  a few  times.  Around  this  disc  is  an  aurora 
of  decreased  brightness  (although  for  large  optical  thicknesses,  it  contains 
most  of  the  energy),  with  width  Bj  on  the  order  of  tens  of  degrees.  Thus, 
increasing  the  f ield-of-view  (FOV)  of  the  optical  receiver  will,  in  general, 
increase  the  number  of  signal  photons  collected.  One  goal  of  these  experi- 
ments is  to  measure  the  signal  increase  with  F0V. 

The  temporal  behavior  of  the  received  signal  with  increasing  FOV  will 
also  provide  needed  information.  The  direct,  unscattered  beam  maintains 
the  same  pulse  shape  as  the  initial  pulse  independent  of  optical  thickness. 
As  the  FOV  of  the  receiver  is  increased,  the  scattered  photons  that  are 
captured  by  the  receiver  may  increase  the  signal  peak  power,  but  will 
probably  also  lead  to  pulse  stretching,  and  to  the  delay  of  the  peak  of  the 
pulse.  Thus,  while  pulse  energy  will  increase,  the  peak  power  may  remain 
the  same,  or  increase  slightly.  The  pulse  car.  be  stretched  to  values  as 
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FIGURE  21.  Mie  scattering  functions  for  various  values  of  the  Mie 
parameter  x=ttD/X.  [24] 
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Fig. 22:  Schematic  Representation  of  Laser  'Angular  Brightness  Distribution 
After  Traversing  a Cloud 


(NOT  TO  SCALE) 


large  as  tens  of  microseconds  [24].  This  will  limit  the  data  rate  to 
tenths  of  megabits  per  second.  More  important,  however,  is  the  fact  that 
the  received  peak  power  is  reduced.  This  factor  is  very  important  to  low 
bit  rate  communications  since  most  conventional  detection  schemes  respond 
to  peak  power  inputs.  In  addition,  the  pulse  stretching  effect  will 
impose  longer  integration  times  in  photon  counting  systems.  For  example, 
if  the  received  pulse  width  increases  from  20  ns  to  20  ys,  the  photon 
counting  interval  must  be  increased  by  10^  in  order  to  collect  all  the 
signal  photons,  and  therefore  increase  the  background  noise  count  also  by 
10  . This  "pulse  stretching"  loss  is  in  addition  to  change  in  background 
noise  resulting  from  increasing  the  f ield-of-view. 

The  transmitter  components  of  the  propagation  link  are  shown  in 
figure  23.  The  laser  is  a General  Photonic  Nd:YAG  laser  emitting  1 MW 
20  nsec  pulses  at  1.06  y with  pulse  rates  to  10  pps . A frequency 
doubling  crystal  converts  about  20  percent  of  this  fundamental  to  first 
harmonic  at  .53  microns.  The  laser  is  triggered  by  a pulse  generator, 
which  will  send  out  a synch  pulse  via  a UHF  transmitter  to  the  receiver. 


The  components  of  the  receiver  are  shown  in  figures  24  and  25.  The 
.53  y radiation  is  collected  by  a variable  FOV  telescope  and  detected  by 
a calibrated  photomultiplier  tube.  The  output  signal  is  digitized  by 
a Biomation  8100  digitizer  and  then  recorded  on  a tape  deck  and  simultane- 
ously displayed  in  real  time  on  a scope.  The  sync  signal  from  the  trans- 
mitter is  detected  and  used  to  arm  the  digitizer  in  preparation  for  the 
incoming  pulse.  The  equipment  for  measuring  the  direct,  nonscattered 
pulse  is  shown  in  figure  25.  The  .53  y radiation  is  collected  by  an  8-inch 
Cassegrain  telescope.  The  telescope  has  a very  narrow  FOV,  and  is  bore- 
sighted  with  the  transmitter  to  detect  only  the  unscattered  beam.  The 
energy  is  again  detected  by  a calibrated  photomultiplier.  A scope  will 
monitor  in  real  time  the  time  history  of  the  radiation.  As  no  pulse 
distortion  is  expected  of  this  direct  beam,  only  the  peak  power  has  to  be 
recorded.  A circuit  to  sense  and  hold  the  peak  voltage  of  the  pulse  will 
then  record  the  value  of  peak  power  of  every  transmitted  pulse.  This  will 
then  allow  determination  of  optical  thickness  of  the  path. 

The  most  important  measurement  in  this  experiment  is  the  time  history 
of  the  pulse  at  the  receiver  as  a function  of  FOV  for  various  optical 
depths  of  fog  conditions.  This  will  yield  the  delay  of  pulse  arrival  time, 
the  amount  and  extent  of  pulse  stretching,  and  the  calibrated  photomultiplier 
tube  will  give  the  flux  density  at  the  receiver,  from  which  the  path  loss 
will  be  calculated  (?/?). 

rec'  trans 

From  the  measurement  of  the  direct,  nonscattered  pulse  power,  an  average 
value  of  single  particle  scattering,  8,  will  be  calculated.  This  value 
will  be  compared  with  that  derived  from  the  particle  spectral  densitv  obtained 
from  the  Knollenberg  counter  and  the  Mie  theory,  and  the  value  that  is 
recorded  from  the  visiometer. 
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Figure  26  is  a schematic  set  of  curves,  normalized  to  clear  air  spreading, 
of  the  data  that  is  expected  from  this  experiment.  They  are  not  to  be  taken 
too  literally,  but  they  should  show  general  trends  and  shapes  that  might  be 
expected.  For  small  values  of  t,  say  t=1,  the  direct  beam  should  dominate,  and 
the  path  loss  should  be  independent  of  FOV.  Consider  the  path  loss  for  t=5. 

For  small  FOV,  the  path  loss  is  e-^,  or  22  dB.  Increasing  the  FOV  will  increase 
the  amount  of  energy  that  arrives  from  the  forward  scattered  component, 

0<FOV>e1,  and  will  probably  saturate  at  FOV  = 6^.  Increasing  the  FOV 
further  will  probably  not  decrease  the  path  loss,  since  for  small  values  of 
the  first  aurora  has  a substantially  larger  value  of  source  brightness 
than  the  outer  aurora.  As  t increases  further,  say  t = 10,  the  brightness 
of  the  first  aurora  decreases,  and  the  multiple  scattering  outer  one 
increases.  There  should  be  an  increase  of  signal  for  values  of  FOV  greater 
than  6^.  Saturation  of  signal  should  occur  at  FOV  angles  in  the  vicinity 
of  02*  For  still  further  increases  of  t,  say  t=  20,  the  position  of  6^  is 
smeared  out  completely.  Increase  of  signal  should  be  observed  smoothly  from 
0 degrees  FOV  to  value  larger  than  FOV  = 30°.  The  main  goal  of  this  experi- 
ment is  the  measurement  of  the  exact  shape  of  the  curves  in  figure  26  and 
determination  of  propagation  regimes. 


As  discussed  above,  the  received  signal  is  also  expected  to  increase 
in  pulse  width  with  increase  of  FOV.  For  0<F0V<6^,  the  increase  is  expected 
to  be  small,  as  the  scattering  angles  are  small.  But  for  0^<FOV<02>  the 
pulse  is  expected  to  stretch  by  orders  of  magnitude.  To  remove  a potential 
ambiguity  in  path  loss  of  figure  26,  we  have  defined  the  path  loss  in  terms 
of  the  peak  power  as  well  as  pulse  energy.  It  will  also  be  very  informative 
to  plot  this  curve  from  the  data,  with  the  vertical  axes  as  the  "Energy 
Path  Loss"  (integrating  the  received  power  over  the  pulse  width  to  obtain 
received  energy).  As  indicated  in  figure  26,  while  increasing  the  field  of 
view  can  substantially  increase  received  signal  energy,  in  pulsed  propagation 
the  increase  of  field  of  view  "may"  only  increase  the  pulse  width  and  not  the 
peak  power,  thereby  yielding  substantially  different  curves  for  the  two  cases. 
The  determination  of  both  these  curves  will  be  important  to  the  system 
designer  of  an  optical  communications  system  operating  through  adverse 
weather . 


3 . 4 Description  of  Atmospheric  Sensors 

Various  atmospheric  remote  sensors  currently  are  in  use  at  or  near 
NELC's  Bldg.  T323  and  are  available  for  assessing  the  propagation  path  for 
the  proposed  measurements.  They  include: 

a.  FM-CW  Radars:  Receives  backscatter  energy  (v  3 GHz)  from  refractive 

turbulent  regions  in  the  lower  atmosphere.  The  humidity  fluctuations  account 
for  essentially  all  of  the  backscatter.  One  radar  is  a fixed-position 
vertically  pointing  unit  located  at  Bldg.  T323  and  the  other  radar  is  a 
mobile  scanning  unit.  These  systems  provide  a means  of  measuring  the 
vertical  extent  of  the  fog  layer. 


b.  Acoustic  Echosounder:  A vertically  pointing  system  (2  kHz)  receiving 

backscatter  energy  from  regions  of  temperature  fluctuations  in  the  lower 
atmosphere.  Using  a transmitter  and  two  receivers,  the  system  is  capable 
of  measuring  horizontal  wind  components  in  a near  real-time  mode  to  an 
altitude  of  500  meters. 
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where  E 


Receiver  Fie Id-of-View  (F.O.V.) 

Pulse  Stretching  Increasing  — 

Figure  26:  Path  loss  versus  f ield-of-view 

tor  various  optical  thicknesses 


c.  Ceilometer  fAN/GMQ  13C)  : A standard  aircraft  runvav  beam  system 

(visible  light,  400  ft.  base  leg)  to  measure  the  height  of  cloud  bases. 

d.  Transmissor.eter  (Ah'/GMQ  10C):  A standard  aircraft  runvav  system 

for  visibility  measurements  using  a visible  beam  of  light  over  a horizontal 
500  ft.  base  leg. 

e.  Visiometer  (MRI  Model  1580A) : A short-path  ('v  2 feet)  forward- 
scatter  instrument  to  measure  visibility  out  to  about  10  km  (with  relative 
humidities  up  to  100%)  . 

f.  Integrating  Nepholometer  (MRI  Model  1550):  An  instrument  for 

measuring  visibility  out  to  about  50  km  (with  relative  humidities  less 
than  70%)  by  measurement  of  light  scattered  through  a sample  of  aerosol 
laden  air. 

g.  Radiosonde  Sets:  One  set  permits  the  measurement  of  winds,  tempera- 

ture and  relative  humidity  from  the  surface  to  above  the  100  mb  level.  The 
other  set  measures  only  temperature  and  relative  humiditv  in  the  same 
altitude  range. 

h.  Standard  Surface  Measurements:  Continuous  records  are  made  of  the 

surface  temperature,  wet-bulb  wind  and  pressure. 


i.  Knollenberg  Drop-Size  Spectrometer:  A system  to  continuously 

measure  and  record  the  drop-size  distribution  of  aerosols  either  on  the 
ground  or  airborne. 

The  location  of  each  of  these  sensors  In  relation  to  Bldg.  T323  is  shown  in 
figure  27.  The  instrument  of  most  importance  to  these  measurements  is  the 
Knollenberg  drop-size  spectrometer.  This  instrument  determines  aerosol-size 
distributions  based  on  theoretical  predictions  of  the  scattering  of  light  from 
a helium-neon  laser  (6328  A)  by  spherical  water  particles.  The  system  outputs 
are  recorded  digitally  on  magnetic  tape  and  are  computer  processed  to  give  the 
particle  number  density  (particles/cm'Vpm  dia.).  Table  I shows  a sample 
computer  print-out  of  a distribution  obtained  during  a dense  Santa  Ana  fog 
which  occurred  on  the  morning  of  25  June  1976. 

3 . 5 Proposed  Schedule 


Dates 


I tem 

Start 

Complete 

1.  Initiate  variable  FOV 
construction 

1 

Oct 

1976 

30  Nov  1976 

2.  Install  neodymium  YAG  laser  at 
Battery  Ashburn 

1 

Oct 

1976 

15  Nov  1976 

3.  Install  data  acquisition  system 
at  Bldg.  T323 

15 

Oct 

1976 

30  Nov  1976 

4.  Initiate  data  acquisition  period 

1 

Dec 

1976 

30  Apr  1977 

5.  Data  analyses  and  final  report 

1 

May 

1977 

30  Sep  1977 
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TABLE  I 

The  measured  density  of  aerosols  at  the  diameters  (urn)  shown 
represent  an  average  of  60  1-second  spectra  taken  over  a 
4-minute  period.  The  visibilitv  during  the  period  as  measured 
by  the  MRI  Model  1580A  visiometer  varied  between  1.18  and 
.179  km.  Also  shown  in  the  table  is  the  forward  extinction 
coefficient  (forward  scattered  plus  absorption)  computed  for 
a wavelength  of  .53  pm  using  the  relation 


s 

r 

N(r)  is  the  measured  distribution  of  spherical  scatterers 
of  radius  r,  Q and  Q are  the  Mie  theorv  absorption  and 
scattering  ef ffcienclls . For  the  example,  the  calculated 


nr  N(r)  0 dr 


(41) 
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•'  is  customarv  in  radiative  t rans 
c ion , so  as  to  consider  the  energv 
solid  angle  d.  and  within  spectral 
making  angle  i with  the  normal  to 
intensity  L,  to  be 


fer  to  make  the  narrow  band  assump- 
dKv  which  flows  in  time  dt  from 
region  ( , v +dv)  across  area  dA 
dA.  I his  defines  the  specific 


I._.  = dE  /dv  dA  cos  $ dfi  dt  . 
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